Entanglement entropy in fermionic Laughlin states 
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We present analytic and numerical calculations on the bipartite entanglement entropy in fractional 
quantum Hall states of the fermionic Laughlin sequence. The partitioning of the system is done 
both by dividing Landau level orbitals and by grouping the fermions themselves. For the case 
of orbital partitioning, our results can be related to spatial partitioning, enabling us to extract 
a topological quantity (the 'total quantum dimension') characterizing the Laughlin states. For 
particle partitioning we prove a very close upper bound for the entanglement entropy of a subset of 
the particles with the rest, and provide an interpretation in terms of exclusion statistics. 
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Introduction — While the renewed attention in quan- 
tum entanglement has mostly been in the context of 
quantum computing, there is now also growing inter- 
est in using entanglement measures for characterizing 
quantum many-particle states. The measure we focus 
on is the entanglement entropy, which is defined by par- 
titioning the system under question into two blocks A 
and B, and using the reduced density matrix of one 
part (e.g., pA = trsp obtained by tracing over B de- 
grees of freedom) to calculate the von Neumann entropy 
Sa — — tr[p^lnp^]. Generally speaking, the entangle- 
ment entropy scales with the size of the boundary be- 
tween the partitions 0. Prefactors, logarithmic cor- 
rections and subleading terms in this basic relationship 
can all contain important physical information concern- 
ing the many-body state. 

By now the entanglement entropy has been studied 
extensively for one-dimensional spin systems, where the 
variation of Sa as a function of the size of block A is a 
sensitive probe into the nature of the ground state 0- 
Entanglement measures have also been shown to con- 
tain revealing behavior in the vicinity of quantum crit- 
ical points [3j and shown to have fundamental implica- 
tions for the efficiency of numerical simulation methods 
for quantum many-body states (Jj. 

Studies of entanglement in higher-dimensional and 
itinerant systems has been less thorough but have now 
also started to attract significant attention (e.g., [ESS 
1^). The case of two dimensions is particularly intriguing 
because of the known existence (and further unconfirmed 
possibilities) of topologically ordered phases. It is natu- 
ral to speculate that entanglement might give a handle 
on topological properties of quantum states, because the 
basic intuition about topologically ordered phases is the 
presence of intricate correlations not easily captured by 
local observables or traditional correlation functions. In- 
deed, Refs. ^] and Q have recently presented the follow- 
ing theorem concerning topologically ordered states: if 
L is the boundary between the two blocks, the entangle- 
ment entropy scales as Sa = oiL — 7 -I- 0{L~^). As usual 
the scaling law applies to situations where A is large and 
the total system is infinite. The quantity 7 (the topo- 



logical entanglement entropy) is the logarithm of a quan- 
tity known as the total quantum dimension. For a state 
with anyonic excitations, the quantum dimensions char- 
acterize the growth rate of the Hilbert space with anyon 
number. For fermionic fractional quantum Hall states in 
the Laughlin sequence with filling fraction v = l/m, the 
topological entanglement entropy is 7 = i Inm. 

In this Letter, we present a detailed study of the en- 
tropy of entanglement in the fermionic Laughlin states. 
These are the best-known many-particle states with topo- 
logical order. We carefully choose useful methods of di- 
viding quantum Hall states into A and B blocks. We 
provide results on the entropy defined by partitioning the 
Landau-level orbitals in the system into blocks, as well 
as that obtained by dividing the particles themselves into 
A and B subsets. The first method allows us, upon ex- 
trapolation to the thermodynamic limit, to extract the 7 
parameter of Refs. [^|^. To our best knowledge, this is 
the first example where the results of Refs. are used 
to compute the total quantum dimension of an experi- 
mentally realized topologically ordered medium, directly 
from a microscopic description. 

Very restricted cases of entanglement in Laughlin 
states, involving one-orbital or one-particle blocks, have 
appeared previously in the literature, in Refs. [3, Il0l| . 
The entanglement entropy in some other topologically 
ordered states has been studied in Ref. iTj . 

Orbital versus particle entanglement — We first ex- 
plain our choice of methods used to partition the Laugh- 
lin system into blocks A and B. The issue is non-trivial in 
the case of itinerant systems where, unlike spin 

systems, the particles are not each fixed in its own site. 
Although this point has not been stressed in the litera- 
ture, conformal field theory results on entanglement scal- 
ing |l2| actually pertain to the blocking of space rather 
than the particles or spins themselves. For quantum Hall 
systems, we find a description of the system in terms of 
the magnetic orbitals more natural than a spatial descrip- 
tion. Thus we will partition the orbitals into two sets A 
and B and calculate the entanglement between them. In 
fact, there is a close relationship between orbital and spa- 
tial partitioning, since the choice of the first Ia orbitals 
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as the A block corresponds approximately to choosing a 
disk-shaped A block with radius proportional to ^/Ta in 
real space. The scaling law of Refs. 0, 1^ thus translates 
to Su^c,VU-l + Oil/lA). 

In addition, we also consider the perhaps more obvi- 
ous method of partitioning the particles themselves. In 
this case, the subsets A and B no longer correspond to 
connected regions in space. We find that the analysis 
of particle entanglement also reveals subtle correlation 
effects. 

Quantum Hall states on a sphere — To describe the 
fractional quantum Hall states, we use the representation 
of Refs. [l^I in which the fermions are placed on a 
sphere containing a magnetic monopole. The magnetic 
orbitals in the lowest Landau level are then represented 
as angular momentum orbitals; for N particles in the 
Laughlin state u = 1/m, the total angular momentum 
is half the number of flux quanta, L — ■^N^ with = 
m{N — 1). The + 1 orbitals are labeled either I = 
to or Lz — —L to +L. The "filling" acquires the 
usual meaning v — N /N^ only in the thermodynamic 
limit. The orbitals are each localized around a "circle of 
latitude" on the sphere, with the / — orbital localized 
near one "pole." For orbital partitioning, we define block 
A to be the first Ia orbitals, extending spatially from one 
pole out to some latitude. In the thermodynamic limit, 
this is equivalent to having a disk-shaped block A] since 
each orbital I is associated with a wavefunction of the 
form z'e~'^' in usual complex coordinate language, a 
disk with I a orbitals has radius ~ ^/Ta- 

In this representation the Laughlin wavefunctions are 
expressed in terms of Schwinger boson operators aj, b\ 
as ni<j("l^] ^ ^1°^])™ |0)- Each term in the expansion 
can be interpreted in terms of orbital occupations. For 
example, in the expansion for the iV = 3, m = 3 wave- 
function, the term involving {aXY {a^)"^ {a\)^ corresponds 
to having the three particles in orbitals 1,2, and 6 respec- 
tively. Together with its permutations, this state reads 
in "orbital" notation [0,1,1,0,0,0,1]. 

Analytic results for orbital partitioning — We first con- 
sider the simplest Laughlin state, m = 3. With only the 
I ~ orbital included in partition A (I a ~ 1), the re- 
duced density matrix pyi is 2 x 2 and diagonal, and has 
eigenvalues N/{3N - 2) and {2N - 2)/(3iV - 2). The 
entanglement entropy converges in the large-TV limit as 
Su=i = ln(3/22/3) + ln22/7iV + OiN-^). 

The case of the first two orbitals in the A block {I a = 2) 
is also exactly solvable. The basis for pA is the set [0,0], 
[1,0], [0,1], [1,1]; however the last one can be dropped 
because the two lowest orbitals are never simultaneously 
occupied in the m — 3 wavefunctions 0. The reduced 
density matrix is still diagonal, with a pair of eigenvalues 
iV/(3A^-2) andasingle (iV- 2)/(3iV- 2). The thermo- 
dynamic limit is Si^=2 = ln3 + A/{9N^) + 0{N-^). 

The case oi Ia = 3 is already more difficult. The re- 
duced density matrix, now of size 5 x 5, is still diagonal; 
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FIG. 1: (Color online.) Entanglement entropies with orbital 
partitioning (5*;^), for various A'^ and m. Inset; maximum 
Sij^ values plotted against square root of positions I'X of the 
maximum, for m = 3 wavefunctions of various sizes. 

the elements are of the form (up to normalization) N, 
a^, UN, {N — QfAr), and {N — 2 — a^)- We find = 1, 
a4 — 92/51 and 015 = 505/203, etc., but a general form 
is not obvious; thus we lack an analytic expression for 
5/^=3 in the thermodynamic limit. 

We now turn to a general Laughlin state of filling 
V = 1/m. For this state there is at most one fermion 
in the first ^{m -\- 1) orbitals ^^; this makes S"/^ ana- 
lytically tractable for I a < ^(to + 1)j as was the case for 
< 2 for the m = 3 case. For Ia < ■^{m + 1), the basis 
state is [0,0,0,. ..,0], [1,0,0,. ..,0], [0,1,0,...,0],... [0,0,0,. ..,1]. 
The reduced density matrix is diagonal. In the thermo- 
dynamic limit the eigenvalues are (m — lA)/m for the 
first state (no particles in A) and 1/m for each of the 
rest. The entanglement entropy is 

Su = -f^)log(^) -/.(^logl) 

which takes the form S*;^ ~ (log m/m)^^! for I a ^ rn. 
Thus for large m the S"/^ versus ^JTa curve starts out 
quadratic, and only later displays the asymptotic linear 
behavior, with a crossover presumably around I a ^ rn. 

Orbital partitioning; numerical results — We now 
present numerical results for entanglement entropy with 
orbital partitioning. The numerical calculation is based 
on Laughlin wavefunctions evaluated exactly using alge- 
braic methods; some details are given later in the article. 

Fig. n shows numerically calculated orbital entangle- 
ments (Si^) for several of our wavefunctions, as a func- 
tion of the number of orbitals {Ia) in the A partition. 
The partition flip symmetry of the entanglement entropy 
(Sa = Sb) implies that S'/^ is arc-shaped with a max- 
imum at l\ = {Nij, + l)/2. The initial increasing parts 
of these curves reflect physics of the macroscopic state, 
while the downward curvature is a finite-size effect. Scal- 
ing information is not immediately obvious because the 
curves bend down before reaching large I a- However, we 
can extract some information about the thermodynamic 
limit by plotting the maximum Si* of the arc-shaped 



3 




^ " , ' No. of fermions, N 

/ ' \ ' \ ' 1 

12 3 4 

1/2 

Sq. root of # orbitals in block, 

FIG. 2: (Color online.) Entanglement entropies in m = 
3 states extrapolated to the thermodynamic limit. Dashed 
line is a fit to —7 + cwJTa giving more weight to the higher 
points. Inset on right plots Sij^ against A'', for I a ~ 1, 3, 6, 10. 
With wavefunctions available upto A'' = 10, the N ^ 00 
extrapolation is reliable for small Ia but is already difficult 
for I A ~ 10 because the Sij^{N) curve is not nearly flat yet. 

curves against the position of the maxima, l\. While 
the detailed finite-size effects are difficult to analyze, it 
is reasonable to assume that the macroscopic functional 
dependence of Si^ (e.g., algebraic or logarithmic depen- 
dence) will also appear in the dependence of the maxi- 
mum. The inset to Fig. ^ shows such a plot for m = 3 
states, with Ia's plotted in a square root scale. The linear 
behavior is manifest. We have observed the same feature 
for m = 5 wavefunctions. We have thus already verified 
the boundary law Si^ oc ^/Ta- 

A more quantitative approach to the thermodynamic 
limit is to do a numerical iV — > 00 extrapolation of the 
Si^ data for each Ia (Fig.^. We fit the Si^iN'^) data 
points (inset to Fig. ^ to functions hke Cq -I- Ci/N°'^ + 
C2/-/V"2, noting that S';^=i_2 indeed have expansions for 
this form. We use various sets {a^} of small integers, 
dropping combinations that give extrapolation functions 
with pathological features at small N~^. The set of val- 
ues obtained by this procedure yields the estimates and 
errors of Fig. 13 

Extracting the topological entropy — We note that the 
—7 + ci^/Ta behavior is expected only for Ia > 3. 

For extracting 7, we thus drop two or more of the lowest 
points from the available S*;^ values. It is not feasible to 
use only the largest values because the extrapolation un- 
certainty is largest there (Fig.EJ. Dropping 2 to 5 of the 
lowest points, and giving various relative weights to the 
higher or lower values, we extract various estimates for 7 
by fitting S*;^ = —7 + ci ^/Ta. The resulting estimate for 
the topological entanglement entropy of the m = 3 state 
is 7 = —0.60 ± 0.13. This is consistent with the theoret- 
ical value of —0.55. However, larger wavefunctions are 
probably required for precisely extracting 7 for models 
where it is not known a priori. 

Particle Entanglement — We now turn to partition- 
ing the Laughlin states into two sets of particles rather 
than two blocks of orbitals. Noting that the entangle- 



ment entropy Sa is largest and equal to the logarithm 
of the number of nonzero eigenvalues when pA has equal 
nonzero eigenvalues, we can set a naive upper bound for 
Sa- When A contains ha particles out of a total of N 
particles in the v = 1/m Laughlin state, the size of the 
density matrix is set in first instance by the number of 
possible ways in which the ua particles can be distributed 
among TV^ -I- 1 = m{N — 1) -I- 1 orbitals; thus 

This bound becomes an exact result for the integer quan- 
tum Hall state m = 1, i.e., the noninteracting case where 
the Laughlin wavefuction reduces to a single Slater deter- 
minant. For the fractional quantum Hall states (m > 3), 
the bound is sharp for ua—I but not for ua > I (Fig.|2J)- 

We note that the absence of terms in the Laughlin 
wavefunctions with adjacent orbital occupancies near the 
poles implies that, for ua > 1, some of the eigenvalues of 
the reduced density matrix are zero, leading to a tighter 
upper bound than above. This idea becomes much more 
effective with the additional observation that the eigen- 
values of Pa are organized in a SU{2) multiplet structure, 
due to Pa commuting with the total angular momentum 
magnitude and z-component (L^^ and i^^) of the ua 
particles in A. When any eigenvalue of pA vanishes due 
to the absence of corresponding states in the Laughlin 
wavefunction, the multiplet structure implies that every 
eigenvalue in the same multiplet must also vanish. 

For UA = 2, m = 3, one observes that the vanishing 
states are part of the L2 — m{N — 1) — 1 = — 1 repre- 
sentation of the SU (2) symmetry algebra. This reduces 
by 2L2 -I- 1 = 2Nif, — 1 the number of 2-particles states 
contributing to the entropy. For ua = 2 and general m, 
there are (m — l)/2 multiplets contributing zero eigen- 
values, leading to a total of \{m — l){2N^ — m + 2) zeros. 
The number of nonzero eigenvalues is 

(^^ + -i("-l)(2iV,-..+2) ^ (^^ + 1 - - 

and the logarithm gives a much better upper bound to 
the UA — 2 entanglement entropies (Fig. I^J. This bound 
is still not sharp because the eigenvalues are not all equal. 
The eigenvalue distribution is itself interesting and may 
be discussed in future work; however, the effect on Sa is 
small. In Fig. Olwe show that the exact values for SnA=2 
rapidly converge to our improved upper bound. 

Generalizing to ua > 2, we have rigorously established 
the following upper bound 

^^<l,g(^A^. + l-(--l)(-A-l)^ . (2) 

We claim that the entropy Sa will be close to this bound 
for UA ^ N. The derivation of Eq. Q relies on a result 
of Ref. for counting degeneracies of Laughlin states 
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m = 5 
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FIG. 3: (Color online.) Entanglement entropy of two par- 
tices (A) with the rest (B). Dots are numerical exact values; 
dashed curves are naive estimates [Eq. Q]; full curves are 
improved bounds taking multiplet structure into account. 



on the sphere in the presence of excess flux, that is, of 
a number of quasi-holes. The combinatorial factor in 
Eq. has the following interpretation: it is the number 
of ways UA particles can be put into Nff, + 1 orbitals, 
such that between two occupied orbitals there are at least 
m — 1 that are empty. Clearly, this counting represents 
a version of Haldane's notion of exclusion statistics [T^ . 

Numerical method — Since it is not known which terms 
in Laughlin wavefunctions are more important for various 
entanglement entropies, we avoided truncation schemes 
and calculated complete wavefunctions. Such calcula- 
tions are in the spirit of combinatorial considerations of 
Refs. rather than that of "exact diagonaliza- 

tion" studies. The enumeration of Laughlin wavefunc- 
tion terms and calculation of coefficients are combinato- 
rially explosive problems [l^l llfij ; symbolic expansion is 
not feasible beyond more that 5 or 6 particles. Instead, 
we made use of results in Ref. where the coefficients 
are expressed as expansions in characters of the symmet- 
ric group, a(^) = E(^) C(A)X(^)- Here (^)'s are integer 
partitions representing terms in the Laughlin wavefunc- 
tion. The conjugacy classes (A) and the coefficients of ex- 
pansion C(^\^ are both found by expanding out a Hankel 
determinant symbolically Characters of the sym- 

metric group Sr ,,,,, , were calculated using publicly 

available group theory code while the Hankel determi- 
nant expansion was done with a symbolic manipulation 
routine coded from scratch. The admissible partitions 
(/i) for Lau ghli n state terms were generated by combina- 
torial rules [la, UM ■ With some effort (50-|- cpu-days for 
the largest wavefunctions), we calculated to = 3 wave- 
functions up to A'^ = 10 fermions and to = 5 wavefunc- 
tions up to TV = 6 fermions. 



Concluding remarks — We conclude that the entangle- 
ment entropy pertaining to the two types of partitioning 
(spatial and particle) reveal different aspects of the intri- 
cate topological order of the Laughlin states. This beau- 
tifully illustrates that, more generally, entanglement en- 
tropy is a valuable probe for quantum order. For the case 
of particle partitioning we remark the following. Trust- 
ing that eq. (2) gives a close bound to Sa, we have for 
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We propose that, in general, the difference between Sa 
and the naive expression (^) based on the number of 1- 
particle states allows a expansion that reveals cor- 
relations and quantum order in a state. 

Our work opens up a number of questions. It would 
be interesting to obtain numeric data for ha > 3 par- 
ticle entanglement (not available at present) to investi- 
gate how well our bound Q works. With further effort, 
some larger wavefunctions may also be obtained, per- 
haps leading to more precise extraction procedures for 
the total quantum dimension. Obviously, other quantum 
Hall states, e.g., those with paired or clustered structure 
and non-abelian statistics, deserve attention, and further 
work in this direction is planned. 

While writing up our results, we learned of parallel 
work |l9l | that has overlap with some issues we consider. 
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